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ABSTRACT 

In this Paper we report on optical spectroscopic observations of the low-mass X-ray 
ly-^ | binary 2S 0921-630 obtained with the Very Large Telescope. We found sinusoidal 

, radial velocity variations of the companion star with a semi-amplitude of 99.1±3.1 

■ km s _1 modulated on a period of 9.006±0.007 days, consistent with the orbital period 

found before for this source, and a systemic velocity of 44.4±2.4 km s _1 . Due to X-ray 
irradiation the centre-of-light measured by the absorption lines from the companion 
' star is likely shifted with respect to the centre-of-mass. We try to correct for this using 

the so-called K-correction. Conservatively applying the maximum correction possible 
O ■ and using the previously measured rotational velocity of the companion star, we find 

I 1 a lower limit to the mass of the compact object in 2S 0921-630 of Mx sin 3 i > 1.90 ± 

0.25 Mq (1 a errors). The inclination in this system is well constrained since partial 
eclipses have been observed in X-rays and optical bands. For inclinations between 
60° < i < 90° we find 1.90 ± 0.25 < M x < 2.9 ± 0.4 M . However, using this 
maximum K-correction we find that the ratio between the mass of the companion 
star and that of the compact object, q, is 1.32 ±0.37 implying super-Eddington mass 
transfer rates; however, evidence for that has not been found in 2S 0921-630. We 
conclude that the compact object in 2S 0921-630 is either a (massive) neutron star or 
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C$ ' a low-mass black hole. 
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1 INTRODUCTION 

The equation of state (EoS) of neutron star matter is in- 
timately related to the physics of the strong interactions 
between fundamental particles and therefore it is of great 
relevance in high-energy and particle physics. Theories on 
the EoS of neutron-star matter at supra-nuclear densities 
provide a firm upper limit on the mass for each EoS. The 
stiffer the EoS, the higher the mass limit. Neutron stars 
with masses well above 1.4 Mq cannot exist for so-called 
soft equations of state, in which matter at high densities is 
relatively compressible (e.g., due to meson condensation or 
a transition between the had ron and quark-gluon phases; 
cf. lLattimer fe PrakashlEoOlh . Therefore, measuring a high 
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mass for even one neutron star would imply the firm re- 
jection of many proposed soft Eo Ss (e.g. see discussion by 
Ivan Paradiis fc McClintocklll995l) . Neutron star masses of 
more than ~ 3 Mq are excluded assuming (among other 
things) that general relativity is the correct theory of grav- 
ity a nd that the velocity of sound i s less than the velocity of 
light llNauenberg fc Chaplineil973t) . Therefore, a lower limit 
on the mass of the compact object of 3 Mq or more implies 
a black hole compact object. 

The most accurately measured neutron star masses we 
have are from precise radio timing measurements of dou- 
ble neutron-star binaries. These empirical mass es cluster 
around 1.35 Mq jThorsett fc Chakrabartvl Il999lh This is 
close to the mass of 1.32 Mq predicted fr om theoretical 
mode l calculations of Type lb supernovae jTimmes et alJ 
1996). However, these double neutron star binaries are 
formed in high-mass X-ray binaries and due to the short 
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life of the massive companion only a limited amount of 
matter can be accreted. In contrast, in low-mass X-ray 
binaries (LMXBs), the rapid rotation and low magnetic 
field of the neutron stars are considere d evidence for much 
larger amounts of accreted matter (see lBhattacharvalll995l 
for a review) . Indeed, from bina ry evolution models, LM XBs 
could accrete up to 0.7 Mq fe.g. lvan den Heuvel fc Bitzarakil 
Il995l) . If so, the masses of their neutron stars could be 
> 2 Mq . Most of the constraints from timing observations 
of the successors of these LMXB primaries, the millisec - 
ond radio pulsars, are weak jThorsett fc Chakrabartvll999l) . 
as the almost purely circular orbits of these systems often 
prevents one from measuring most of the post-Newtonian 
parameters that underlie the precise masses for the 
double-neutron star binaries. With masses of 1.57±°;f Mq 
and 2.2±0.4M Q , the msec pulsars PS R B1855+09 and 
PSR J0751+180 7 are exceptions to this (iNice et al.ll2003t 
iNice et al.ll200l. 95% conf.). 

Dynamical mass estimates or dynamically determined 
limits on the mass of the neutron star in LMXBs are avail- 
able for the two transient sources Cen X-4 and XTE J2123- 
058, the Z-source Cyg X-2, and t he pulsar 2S 1822- 
371 (Cen X-4; M NS = 1.3±0.6 Mq IShahbaz et al1ll993t 
XTE J2123-058; M NS = 1.55±0.31 ~M^ ICasares et all 
| 2002l: Cve X-2: Mns = 1.78±0.23 M^ ICasares et alJll998t 
lOrosz fc Kuulkerj Il999|: 2S 1822-371: Mm > 1.14±0.06 
Mr^ lJonker et all2003tlcasares et al.l2003h . The problem is 
that many of the companion stars are small, late type stars 
which at a distance of more than a few kpc are too faint 
to observe, especially since the integration times are limited 
to typically l/20th of the orbital period in order to avoid 
Doppler smearing of the spectral lines due to the binary 
motion. Furthermore, the binary inclination is often poorly 
constrained and effects of irradiation shift the centre-of-light 
with respect to the centre-of-mass. 

Some of the aforementioned problems will be less se- 
vere when observing long orbital period systems for which 
the inclination is well constrained from the fact that eclipses 
are observed: 2S 0921-630/V395 Car is such a high in- 
clination LMXB with a long orbital period (P or j, ~ 9.0 
days; e.g. JCowlev^3^dJ^)82|, iBranduardi-Ravmont et all 
119831. and iMason et alJll987f) . The inclination of 2S 0921- 
630 is relatively well constrained and must be high 
since partial eclipses of the compact object and accre- 
tion disc have been o bserved in both X-rays and op- 
tical wave bands (e.g. [Braiidjiardi^^mnont et all Il98lt 
IChevalier fc Ilovaiskvll98ll : lMason et alll987Tl . 2S 0921-630 
is thought to be a halo object in orbit with a K0-1 III com- 
panion st ar; absorption lines of the companion have been 
detec ted iBranduardi-Ravmont et all Il983l : IShahbaz et all 
1999). The corresponding radial velocity curve would pro- 
vide a limit to the mass of the companion star. In this 
Paper we report on spectroscopic observations of 2S 0921- 
630/V395 Car. 



2 OBSERVATIONS AND ANALYSIS 

We observed 2S 0921-630/V395 Car with the FORS2 spec- 
trograph mounted on the Yepun Very Large Telescope 
(VLT). In the period Dec 26, 2003 - Mar 13, 2004, 22 spec- 
tra using the 1200R+93 and 22 spectra using the 1028z+29 



grating have been obtained in Service Mode. The exposure 
time was 1300 s for each observation and a slit-width of 
0.4" was used on each occasion. The dispersion was 0.75 
A per pixel for the spectra obtained with the 1200R+93 
grating and 0.86 A per pixel for the spectra obtained with 
the 1028z+29 grating. With a slit-width of 0.4" the two- 
pixel resolution is approximately 65 km s _1 at 6500 A. After 
each observing night Helium-Neon-Argon lamp wavelength 
calibration spectra were obtained. The seeing varied between 
0.5" and 1.4" and the spectra were obtained when the source 
had an airmass of ~ 1.2 — 1.4, except on one occasion on Jan 
20, 2004 when the airmass was nearly 2. 

The spectra were bias subtracted, Hatfield corrected, op- 
timally extracted and wavelength calibrated using the iraf 1 
reduction package. The rms scatter in the wavelength cali- 
bration was «0.03A. We further applied a small correction 
to the wavelength solution by aligning sky lines observed in 
each of the spectra. Due to variations in seeing, resulting in 
variable slit losses, the signal-to-noise ratio varied between 
the spectra from «80-240 over the wavelength range 5920- 
6520 A. Next, we exported the extracted spectra to the data 
analysis package MOLLY. 

In molly the observation times were first corrected to 
the Heliocentric Julian Date time frame (using UTC times) . 
Next, we normalised and rebinned the spectra to a uniform 
velocity scale removing the Earth's velocity. Below, we first 
discuss results obtained using the 1200R+93 grating spec- 
tra. We cross-correlated the 1200R+93 grating spectra with 
spectra of template stars rebinned to the same velocity scale. 
Nine template star spectra with spectral types ranging from 
G5-K7 had been obtained with the Keck telescope at similar 
resolution. We fitted a sinusoid to the measured velocities 
as a function of time. The fit-parameters are: the phase, the 
semi-amplitude, and the period of the sinusoid and the sys- 
temic velocity. However, since the errors on the velocity rep- 
resent the statistical error bars only, they do not include the 
systematic effects which we found to be important (we will 
come back to this in the Discussion). Therefore, the reduced 
X 2 of the fit was much larger than 1. In order to estimate 
realistic errors on the fit parameters, we increased the size 
of the error bars such that the reduced \ 2 was close to 1. 
We found a best-fit period of 9.006±0.007 days, consistent 
with the orbital period of «9 days found before (here and 
below we give 1 a error bars). From this we derive an orbital 
ephemeris of HJD = 2453000.49(8) + Nx9 d .006(7) (UTC) 
where phase zero is defined as superior conjunction of the 
neutron star and the number in between brackets denotes 
the uncertainty in the last digit. The derived ephemeris does 
not significantly depend on the spectral type of the used 
template star and i s consistent within the errors with the 
ephemeris given by IMason et all il987T) . We determined a 
radial velocity semi-amplitude of 99.1±3.1 km s _1 using the 
K1V template star HD 124106 (see Figure[TJ. The radial ve- 
locity semi-amplitude does not significantly depend on the 
spectral type of the template star. I.e. the amplitudes we 
found varied by 2 km s _1 between the minimum and maxi- 
mum amplitude for template spectral types in the range of 
G5-K7. We found a systemic velocity of 44.4±2.4 km s . 
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Figure 1. Radial velocity of the companion star in 2S 0921-630. 
Two cycles have been plotted for clarity. Phase zero is superior 
conjunction of the neutron star. Overplotted is the best-fit sinu- 
soidal radial velocity curve. The data points are shown with their 
formal statistical errors only. 

After correcting each individual spectrum for the ob- 
served sinusoidal velocity shift appropriate for the orbital 
phase we created an average spectrum in the frame of the 
donor star for the data of 2S 0921-630. One can in principle 
derive the rotational broadening from an opti mal subtrac- 
tion (see for the description of this technique iMarsh et alJ 
1994). However, our data is not well suited for this since the 
dispersion at which both the template star and the source 
were observed is relatively low (~ 30 km s _1 per pixel in 
the 5920-6520 A spectral range for the template star, this 
spectrum was resampled to the dispersion of ~ 34 km s~ 
per pixel in the 5920-6520 A spectral range at which the 
source was observed) . Since the template stars were not ob- 
tained with the same telescope and instrumental set-up, the 
instrumental profile will be different from that of the object 
spectra; together with the low resolution this means that it 
is not possible to reliably measure the rotational broaden- 
ing when it is of the order of the spectral resolution. Earlier 
results o n the rotational broa dening of 64±9 km s _1 mea- 
sured byJShahba^t^JjEH^) showed that this is likely the 
case. IShahbaz et alJ (Il99<t) found a best-fit spectral type of 
K0 for the companion star; we found that the residuals were 
smallest when we used the Kl template star HD 124106 
(irrespective of broadening). We further found that on av- 
erage 20 per cent of the light comes from the companion 
star in the 5920-6520 A range. In Figure [5] (Left panel) we 
plot the average spectrum in this range, the broadened tem- 
plate star veiled by a source of constant light contributing 80 
per cent of the light, and the residuals after subtraction the 
veiled template star spectrum. We used the spectral range 
from 5920-6520 A for the optimal subtraction, however, the 
range in wavelengths between 6250-6350 A was masked since 
it contains a strong interstellar absorption feature. 

We further examined whether the contribution of the 
Kl companion star to the total amount of light in the 5920- 
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Figure 3. The fractional contribution of the Kl companion star 
to the total light in the 5920-6520 A wavelength range as a func- 
tion of orbital phase. Two cycles have been plotted for clarity. 
Phase zero is superior conjunction of the neutron star. 



6520 A range varies as a function of the binary orbital phase. 
We averaged spectra obtained at the same phase using phase 
bins of width 0.1. We optimally subtracted the template 
spectrum of the Kl star from the average spectra in each 
phase bin. The measured fractional contribution of the com- 
panion star to the total amount of light in the 5920-6520 
A range is plotted in Figure |3 The large difference in frac- 
tional contribution from the K-star between phase and 
phase 0.1 is difficult to explain. 

Besides the stellar absorption lines apparent in the 
wavelength range 5920-6520 A shown in the top spectrum of 
the Bottom left panel in Figure [5] additional stellar absorp- 
tion lines were present in the range 6520-6800 A (see the 
Top left panel in Figure Furthermore, several emission 
lines are present in the spectrum. Unfortunately, our obser- 
vations lack the resolution to resolve these accretion disc 
lines and he nce only one broad em ission peak was observed 
whereas e.g. lShahbaz et all il999i) observed that the Ha and 
He I emission lines at 6562.76 A and 6678.15 A, respectively 
have the double peaked profile typically observed in accre- 
tion disc spectra. We detected a broad emission peak at 
^6820 A which we could not positively identify. We checked 
whether CCD defects were present at the position relevant 
for this wavelength in the flatfield images or in the individual 
images but none were found. 

So far, we have only discussed spectra obtained with the 
1200R+93 grating. In Figure|2] (Bottom right panel) we show 
the phase folded average spectra of the 1028z+29 grating in 
the range between 8300-8900 A. Clear donor star features 
such as the Ca II triplet at 8498 A, 8542 A, and 8662 A, 
and the Mg I line at 8806.75 A can be seen at orbital phases 
where the back, non-irradiated side, of the companion star 
is observed. Furthermore, superimposed on the Paschen disc 
emission lines, Paschen absorption lines can be seen. Unfor- 
tunately, the combined effect of two unresolved broad ac- 
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Figure 2. Top left panel: The spectrum after averaging all the spectra obtained with the 1200R+93 grating. Top right panel: Same as 
the top left panel but using spectra obtained with the 1028z+29 grating. Bottom left panel: From top to bottom: The average spectrum 
of 2S 0921-630 in the range 5920-6520 Ain the frame of the companion star. The Kl template star (HD 124106) spectrum veiled by an 
accretion disc contribution of 80 per cent. The residual spectrum after subtraction of the veiled template star spectrum. For clarity, the 
spectra have been normalised and shifted vertically. Bottom right panel: Average spectra as a function of phase, the number on the right 
indicates the centre of the orbital phase bin (there are no spectra in the phase bin from 0.35-0.45 obtained with the 1028z+29 grating 
hence only 9 spectra are plotted). The Call and Mgl lines from the companion star have been identified below the average spectrum at 
phase 0.0. 



cretion disc emission lines, a Paschen absorption line and, 
at some wavelengths, and at certain orbital phases, a Ca II 
absorption line is too difficult to disentangle with the lim- 
ited resolution data in hand. Therefore, we did not use these 
spectra for our radial velocity study. Besides the features vis- 
ible in the bottom right panel of Figure [5] we identified two 
other lines in the 1028z+29 grating spectra; a weak Fe I ab- 
sorption line complex near 7750 A and a Paschen emission 
line at 9229 A(see the top right panel in Figure 



3 DISCUSSION 

We determined the radial velocity curve of the companion 
star in 2S 0921-630/V395 Car using VLT/FORS2 spectro- 
scopic observations. A sinusoidal fit to the radial velocity 
measurements with a semi-amplitude, K2, of 99.1±3.1 km 



s _1 and with a systemic velocity, 7 =44.4±2.4 km s _1 rep- 
resents the data well. In these error estimates we artificially 
increased the error bars on the individual velocity measure- 
ments by a factor 4 in order to obtain a reduced v 2 of the 
fit close to 1. Before we discuss the implications of the K- 
velocity, we investigate possible reasons for the large reduced 
X 2 obtained when using the formal statistical errors only. 

Near quadrature at phase 0.25 the velocity measure- 
ments obtained on different nights differ well beyond the 
statistical errors. From Figure [3] we can see that the com- 
panion star contribution to the total amount of light in the 
5920-6520 A range varies as a function of orbital phase. The 
phasing of this variation suggests that the centre of light 
emitted by the K-star is likely to be shifted with respect 
to its centre of mass. An explanation is that the inner side 
of the companion star is heated by the (X-ray) irradiation 
coming from (near) the compact object, reducing the equiv- 
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alent width of t he K-star stellar lines in the range 5920-6520 
A (see also e.g. IWade fc Hornelll98l) . This manifests itself 
in several ways: 

(i) a lower fractional contribution to the total amount of 
light from the K star at e.g. phase 0.5 when using the opti- 
mal subtraction technique (e.g. Figure 

(ii) an increase in velocity (especially near quadrature) with 
respect to the velocity associated with the centre of mass of 
the companion star, since the line profiles will, as a result 
of the stellar rotation, be preferentially shifted to either a 
larger red- (phase 0.25) or blue-shift (phase 0.75). Hence, 
variations in the amount of (X-ray) irradiation between the 
observations could explain the variations in velocities near 
quadrature as observed in Figure Since the observations 
were obtained with at least one day in between, with a base- 
line of several months, it seems likely that the irradiation 
changed between the observations. 

Furthermore, the velocity measurements around phase 
0.5 do not fall on the best-fit sinusoid. Those at phases 0.4- 
0.5 lay above the best-fit sinusoid, whereas those at phases 
0.5-0.6 fall below the cur ve. This could well be explained 
by the Rossiter effect (cf. iHilditchl 1200 ill . Such deviations 
occur if the companion star is partially eclipsed by the ac- 
cretion disc. Such a partial eclipse of the companion star 
can also help explain the low fraction of light coming from 
the companion star at phase 0.5 (see Figure 0. The radial 
velocity measurements at phases just prior to phase 0.5, the 
partial eclipse of the companion star, are biased to higher, 
red shifted velocities since the part of the companion star 
that is eclipsed at those phases is mainly rotating towards 
the observer. Hence, the rotational broadening of the line 
will be biased towards the red. The reverse holds for phases 
just after the partial eclipse. 

For the rotational broadening of the stellar lines we tak e 
v sin i = 64 ± 9 km s^ 1 as found bv lShahbaz et alJ ((1999). 
For a Roche lobe filling companion star the following re- 
lation holds: ^ = 0.46[(1 + qfq]^ fe.g. IWade fc Hornd 
Il988l : q is denned here as M2 divided by the mass of 
the compact object, Mx). Before we can determine q, 
and hence with K2 the mass of the compact object, we 
have to take effects of the non-uniform absorption distri- 
bution on the stell ar surface caused by ir radiation into ac- 
count. To this end, Wad e fc Hornd lll988T) developed a so- 
called K-correction. Following the procedure laid out in 
IWade fc Hornl (|l988h . K 2 , CO rr = K obs - AK. K 2 ,corr is the 
corrected observed (K t s ) radial velocity semi-amplitude 
and AK — fRiKija-i, = fvsini. Here a 2 is the distance 
between the centre-of-mass of the binary and the centre-of- 
mass of the companion star and / is a geometrical correc- 
tion factor less than 1. In the extreme case that the hemi- 
sphere facing the compact object does not contribute to 
the observed stellar absorption lines at all and the other 
hemisphere has a uniform absorption line strength, f = 
JWade fc Horndll988h . Since we still have a contribution of 
the Kl star of around 7 per cent at phase 0.5 we know that 
this is an overestimation of the K-correction in the case 
of 2S 0921-630. Hence, applying this upper limit on the K- 
correction, AK ^ 26.9±3.8 km s , we get a firm lower limit 
of K2, C orr 72.2 ± 4.9 km s _ . Plugging-in the numbers 
given above yields q=1.32±0.37; n ote that in derivi ng this 
we used the equation for Rija of iPaczvnskil il97ll) which 
is valid for < q < 0.8. However, the difference with the 



equation of lEggletonl dl983l) is small for q<1.7 and does not 
change the result significantly given the large uncertainty in 
v sin i. The uncertainty in q is dominated by the large error 
on the rotational velocity. Thus, we find for the minimum 
mass of the compact object, Mx sin i > 1.90 ±0.25 Mq. As 
mentioned in the Introduction, the inclination of 2S 0921- 
630 must be high. Furthermore, we have presented evi- 
dence for a partial eclipse of t he companion star. For in- 
clinations, 60° < i < 9 0° (cf. IfBranduardi-Ravmont et alJ 
119831: iMason et alJll987l) the mass of the compact object is 
1.90 ±0.25 < M x < 2.9±O.4M . 

Using vsini — 2nR2 sini/P, we get for the radius of 
the companion star, 7?2, -R2sini = 11.4 ± 1.6 _Rq. Under 
the assumption that the companion star fills its Roche lobe, 
the mass of the donor star, M2, depends on ly on its ra- 
dius and the orbital period fcf. |Paczvnskill9'7lh . which gives 
M 2 sin 3 i — ^234 3 x — (h ere trie orbital period is measured 
in hours, and the mass and radius of the companion star 
are in solar units). This gives M2 sin 3 i — 2.5 ± 1.1 Mq. 
This M2 suggests an interesting evolutionary stage of the 
binary. M2 is more massive than standard donor stars in 
LMXBs. Recently, t here have been a number of evolution- 
ary c alculations fe.g. lTauris et alfeOOftlpodsiadlowski et alJ 
2002) showing that systems with donor masses as high as 
~5 Mq can enter a relatively long phase of mass transfer 
to a neutron star after the donor has lost most of its mass. 
2S 0921-630 would nicely fit in this scenario, evolving along 
the track of an initial 3-4 Mq donor. However, in this evo- 
lutionary scenario the current time-averaged mass-transfer 
rate is expected to be high, certainly super-Eddington. This 
is at odds with the finding from studying X-ray lines as 
measured with Chandra and XMM-Newton that the ob- 
served X-ray luminosity of ~ 10 36 erg s _1 (for a distance 
of ~7 kpc) is likely clo se to the intrinsic source luminos- 
ity llKallman et all2003r) . Perhaps the current mass-transfer 
rate is much lower than the time-averaged transfer rate, or 
the donor mass and/or the mass-ratio are overestimated. 

Above, we assumed that the companion star is in co- 
rotation with the binary orbit. Synchronization of the ro- 
tational period and the binary period occurs on timescales 
shorter than the circularization timescale fe.g. lTassoull988T) 
and no signs for orbital eccentricity have been found (see 
Figure 0. Furthermore, using the c ompanion star parame- 
ters derived above and equation 7 of lTassoul <ll988l) we find 
that the synchronization timescale is <J 100 yr for 2S 0921- 
630, justifying this assumption. We conclude that the com- 
pact object in 2S 0921-630 is likely to be a massive neutron 
star or a low-mass black hole. If indeed the compact object 
in 2S 0921-630 is a massive neutron star the mass measure- 
ment would rule out soft equations of state. Future high 
resolution optical observations should provide a much more 
accurate handle on the rotational velocity of the companion 
star, yielding a more accurate constraint on the mass of the 
compact object. 
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